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1. INTRODUCTION 

Due to the damage caused to our environment as a result of the widespread use of petroleum 
factories and cars, attention is turning to the use of renewable and environmentally friendly energies such as 
solar energy. This latter depends on PV panels to convert sunlight into electricity. The market has many types 
of PV panels. The difference lies in the material from which it is made and the size, in addition to the number 
of PV cells used. Each PV panel has its own data sheet, which does not include the values of the photocurrent 
source, the resistances, nor the quality factors of the diode and its reverse saturation current. These 
parameters are required to adjust the PV panel models. Three equivalent circuits can be used for the PV 
panel, namely: the one-diode model, the double-diode model, and the three-diode model. Our work is 
focused on the two-diode model. Many techniques have been used and developed for the two-diode model of 
solar cell/module during the last few years. Among these techniques are: the salp swarm algorithm [1], the 
biogeography-based heterogeneous cuckoo search algorithm [2], the performance-guided JAYA algorithm 
[3], the hybridized interior search algorithm [4], the sine cosine algorithm [5], the flexible particle swarm 
optimization algorithm [6], a modified flower algorithm [7], and the Lozi map algorithm [8]. The particle 
swarm optimization algorithm based on anarchic universal heterogeneous learning [9], the cuckoo search 
algorithm and grey wolf optimizer [10], the interval branch and bound algorithm [11], a multitask 
optimization technique that is evolutionary [12], the levy flight backtracking search algorithm [13], the 
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Nelder-Mead moth-based orthogonal flame technique [14], the farmland fertility optimization algorithm [15], 
and the differential vectors reusing backtracking search algorithm [16]. 

The particle swarm optimization algorithm based on multi-swarm spiral leader [17], the algorithm for 
Laplacian Barnacles mating based on neighborhood strategy [18], a logistic chaotic JAYA algorithm [19], the 
slime Mould algorithm [20], a differential evolution algorithm based on an ensemble of self-adaptive [21], the 
salp swarm algorithm based on opposition learning [22], the improved equilibrium optimizer algorithm [23], 
and the Harris Hawks optimization algorithm [24]. The adaptive teaching-learning algorithm [25], the particle 
swarm algorithm based on parallel computing of the niches [26], the vertical and horizontal crossover of the 
Nelder-Mead Simplex and Harris Hawks algorithms [27], the ant lion optimizer algorithm [28], the particle 
swarm optimization technique based on the ensemble of the fractional chaotic [29], the backtracking search 
algorithm [30], the JAYA algorithm [31], and an optimization algorithm for generalized normal distribution 
[32]. The equilibrium optimizer algorithm [33], the distribution algorithm [34], a shuffled frog-leaping 
algorithm [35], the reinforcement learning differential evolution algorithm [36], the spherical evolution 
algorithm [37], the supply-demand-based optimization algorithm [38], an improved slime Mould algorithm 
[39], and the combined self-adapting algorithm with Newton-Raphson method [40]. The non-uniform 
mutation genetic algorithm [41], the marine predators’ algorithm [42], the frog-leaping algorithm [43], the 
Rao-1 algorithm [44] and its logistic chaotic optimization algorithm [45], an advanced slime Mould algorithm 
[46], an enhanced algorithm to optimize the adaptation of the butterfly [47], and a combination of two 
algorithms: experience learning with teaching-learning [48]. The whippy Harris Hawks Optimization 
Algorithm [49], a combination of two algorithms: the adaptive grey wolf optimizer and its chaotic algorithm 
[50], a combined shuffled frog leaping algorithm with a driven ensemble by multi-strategy [51], an enhanced 
gaining-sharing knowledge algorithm [52], an enhanced equilibrium optimization algorithm [53], the Laplace 
cross search algorithm has been combined with the Nelder-Mead simplex algorithm [54], the improved JAYA 
algorithm [55], and a combination of two algorithms: gradient-based optimization and random learning [56]. 
In addition to the arithmetic optimization algorithm [57], the heap-based optimizer algorithm [58], the 
heterogeneous differential evolution algorithm [59], the bio-inspired metaheuristic algorithms [60], and a 
combination of two algorithms: Crisscross and Nelder-mead [61]. In this paper, the MATLAB/Simulink 
simulator has been used to build the proposed mathematical equations for the identification of the parameters 
of the solar panel. Per-unit representations of the current-voltage and power-voltage characteristics of a PV 
panel have been represented at and NOCT. Moreover, unknown parameters have been calculated for three 
atmospheric conditions, namely: STC, NOCT, and random atmospheric conditions. 


2. PARAMETER IDENTIFICATION OF PHOTOVOLTAIC PANEL 

The basic components of a solar panel are the solar cells. Therefore, the modeling of a solar panel is 
based on the solar cell’s equivalent electrical circuit. Figure 1 represents the equivalent electrical circuit of 
the two-diode model for a solar panel. This circuit consists of a photocurrent source, [,,,p, in parallel with the 
two diodes, D1 and D2, in addition to the parallel and series resistances, Rspp and Rsp, respectively. The 
characteristic equation of the PV panel’s current has been obtained using Kirchhoff's laws with consideration 
of Shockley’s diode equation, as given by (1). This equation has seven unknown parameters: the 
photocurrent source, Ippp, the two diode reverse saturation currents, [gp,and Igp2, the two thermal voltages, 
Vrnp1» Vrnp2. the series resistance, Rsp, and the parallel resistance, Ro,p. In this paper, MATLAB/Simulink 
software has been used to implement the proposed mathematical equations for simultaneously calculating 
these seven unknown parameters at different atmospheric conditions. The thermal voltages can be determined 
using (2) and (3). np; and np2 are the quality factors of the two diodes. Ky is the constant of Boltzmann 
(1.38065*10°J/K), q is the magnitude of the electronic charge (1.602*10°!°C), T, is the ambient temperature 
(K), Nsis the number of solar cells connected in series, and Tstc is the temperature at STC (298K). 


Rsplp +Vp Rsplp +Vp icetaate 

Ip = Ipnp — Topi (exp VTnP1 — 1) — Ippo (ex VrnP2. — 1) — Pa sprP ; spp (1) 
ShP 

Vrnpr = “**T, Ns (22) 0.7 (2) 

Vrnp2 = “2*2T, Ns (2) 0.7 (3) 


The short-circuit current and the open-circuit voltage can be calculated using (4) and (5), respectively. 
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The temperature coefficients of the short-circuit current and the open-circuit voltage, ki and ky, 
respectively, are usually given in the data sheet of a PV panel. JR is the solar irradiation at an ambient 
temperature, and /Rsrc is the solar irradiation at STC (1000 W/m7). k;;¢ has to be adjusted to obtain the value 
of the short circuit current. Moreover, the value of the open circuit voltage can be obtained by varying Kyoc. 
The two saturation currents’ values can be calculated using (6) and (7). Moreover, the value of the 
photocurrent source can be obtained using (8). 
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The series and parallel resistances are given by (9) and (10), respectively. The number of the parallel 
solar cells is Np. S is the surface of the solar panel. The resistance rs must be varied until the maximum 
power point (MPP) of the solar panel is attained at STC. For any other atmospheric condition, kgs has to be 
varied to attain the PV panel’s MPP. 


P 
Bie ee eg (5max) (“2 P ae (9) 
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Ns (5 ) Ta (- 72) 
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Figure 1. Equivalent electrical circuit of the double-diode model for PV panel 
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3. RESULTS AND DISCUSSION 

In this section, the MATLAB/Simulink simulator has been used to implement (1)—(10). It has the 
ambient temperature, solar irradiation, and the data sheet as inputs, while its outputs are the current and 
power of the PV panel. Different sizes of PV panels for different markets have been used, such as Suntech of 
SOW, SW 85 poly RNA of 85 W, and JAM60S01 of 320 W. The parameters of these PV panels have been 
obtained at STC and NOCT. Moreover, the parameters of the Suntech PV panel have been obtained under 
random atmospheric conditions. 


3.1. Results at STC and NOCT 

We achieved the MPP of each used panel by varying only the value of rs for the series resistance 
given by (9) at STC. We adjusted the variable kp,of the series resistance for each used PV panel to achieve 
its maximum power point under NOCT (Gr = 800 W/m’, Ta = 20 °C). The data sheets of the used panels are 
given in Table 1. Table 2 shows the adjusted values of rs and kp,in addition to the calculated parameters. 
Figure 2 shows the per-unit curves of I-V and P-V for the used panels (a) at STC and (b) under NOCT. Pds 
and Ids are used for the per-unit values of the power and the current, which are given in the data sheet, 
respectively. There are three highlighted points for each used panel. The per-unit points of the short circuit 
current and open circuit voltage are equal to one, so their calculated values are equal to the data sheet's 
values. The third point is the maximum power point. Its per-unit values are very close to one. Therefore, the 
calculated values are very close to the data sheet values. This closeness is calculated using (11). Table 3 
shows the obtained values for the absolute relative error. 


|calculated value—data sheet value| 
Ae % = « 100 (11) 


data sheet value 


Table 1. Values from the data sheet of PV panels 
STC conditions NOCT conditions 
SOW 85 W 320 W 50 W 85 W 320 W 
PratW) = 50.75 85.204 = 320.0688 + =36.192 62.103 235.162 


Tpmax (A) 2.9 4.76 9.56 2.32 3.81 7.66 
Vemax(V) 17.5 17.9 33.48 15.6 16.3 30.70 
Tsc(A) 3.2 5.02 10.5 2.59 4.04 8.05 
Voc (V) 21.8 22.1 40.8 19.8 20.1 37.61 


Table 2. Calculated parameters of PV panels 


STC conditions NOCT conditions 

50 W 85 W 320 W 50 W 85 W 320 W 
Kise / / / 0.92687 0.9571 0.98245 
Kvoc / / / 0.5875 0.5858 0.495585 
Krs / / / -2.8505 0.17839 0.1520 
r(Q) 6.645 10.0085 9.6265 6.645 10.0085 9.6265 
Vinmi(V) 1.6827176 0.8737 1.4022 1.3235872 0.68724723 1.102989 
Vinm2(V) 1.6827176 0.8737 1.4022 1.3235872 0.68724723 1.102989 
Isc (A) 3.2 5.02 10.05 2.59000257 4.04006589 8.0500727 
Iph (A) 3.20016593 5.0220 10.0540 2.590746323 4.04249936 8.0547659 
Voc(V) 21.8 22.1 40.8 19.800792 20.10003 37.6100062 
Rs (mQ) 92.40804 323.22097 309.8726 515.244988 466.28187 452.784478 
Rsh(Q) 1797.7331 774.12 776.6496 1797.73315 774.1221 776.6496 


Thani 10°A) 3.76797 25.84* 10° — 1.1558*10°% 0.4106987 3.9896936* 107  6.2190.3391* 10° 
Tuarx(10°°A) 3.76797 25.84*10% — 1.1558*10°% 0.4106987 3.9896936*107  6.2190.3391*10° 
Témax(A) 2.91 4.748 9.5669 2.3604 3.8334 7.7003 
Vemax(V) 17.44 17.9452 33.456 15.3329 16.2006 30.5393 
Puimax(W) 50.7504 85.2038096 320.0702 36.19177716 62.10338004 235.16177179 


Table 3. Absolute relative errors of PV panels 


STC conditions NOCT conditions 
50 W 85 W 320 W 50 W 85 W 320 W 
Aej5¢ % 0 0 0 0.0000995 0.00163094 0.0009031 
Aeyoc % 0 0 0 0.004 0.00014925 0.00001648 
A€imax % 0.34482 0.2521 0.072175 1.741379 0.61417 0.52610966 
Aeymax % 0.34285 0.2525 0.071684 = 1.722179 0.6131288 0.52345 


So 


Aepmax % 0.000788 0.000223 0.0004 0.0006157 _0.00061195 _0.00009704374 
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3.2. Results at arbitrary atmospheric conditions 
One experiment had been done at the university using a Suntech PV panel. The atmospheric 
temperature was 17 °C and the solar irradiation was 1176 W/m?. By taking atmospheric circumstances into 
account, we adjusted the variables, Kp skvoc and kisc, to make the characteristic curves as close as possible to 
the measured curves. Therefore, the parameters are obtained simultaneously. Figure 3 shows the 
characteristic curves of the Suntech panel at arbitrary atmospheric conditions. The adjusted values of the 
variables and calculated parameters are given in Table 4. The measured and calculated data for the three 
points, in addition to the resulting absolute relative error, are given in Table 5. 
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Figure 2. Per-unit curves of I-V and P-V for the used PV panels: (a) at STC and (b) under NOCT 


Figure 3. Suntech solar panel: I-V and P-V characteristics at arbitrary atmospheric conditions 


Table 4. Calculated parameters at arbitrary atmospheric conditions 


The parameter Calculated value The parameter Calculated value 


Kise -0.22 Vinmi (V) 1.1135299 
Kvoc -0.95 Vinmz2 (V) 1.1135299 
Rrs -6.7444 Teas (DA) 32.7536 
Rs (mQ) 521.87354 Tan (A) 32.7536 
Rsh (Q) 3109.020864 Iph (A) 3.0651783 
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Table 5. PV panel under random atmospheric conditions 


Isc(A) Voc (V) Tpmax (A) Vmax (V) Pinax(W) 
Measured data - - 2.811 15.408 43.311888 
Calculated data 3.0646637 19.663987 2.8321 15.2933 = 43.31215493 
Ae % - - 0.7506 0.7444 0.000616 


4. CONCLUSION 

The goal of this paper is the calculation of the parameters for the photovoltaic panel with a two- 
diode model. Therefore, mathematical equations have been proposed and implemented using the 
MATLAB/Simulink simulator. Per-unit representations of the current and power of the solar panel with 
comparisons were obtained for STC as well as NOCT. Moreover, comparisons have been made under 
arbitrary meteorological circumstances. For the three cases, the absolute relative error of the MPP did not 
reach 0.0008%. 
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